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ABSTRACT 

Both the power spectrum of the Cosmic Microwave Background and the Big Bang Nucleosynthesis 
theory combined with observations of light elements imply that the baryon density is about 4.5% of 
critical. Most of these baryons are observed at redshifts greater than two, but have remained elusive at 
lower redshifts. Hydrodynamic simulations predict that the low redshift baryons are predominately in 
a warm- hot intergalactic medium (WHIM) , which should exhibit absorption and emission lines in the 
soft X-ray region. The WHIM is predicted to be most dense around clusters of galaxies. We present 
our XMM-Newton RGS observations of X Comae, an AGN behind the Coma cluster. We detect 
absorption by Ne IX and O VIII at the redshift of Coma with an equivalent width of 3.3 ± 1.8 eV and 
1.7 ± 1.3 eV respectively (90% confidence errors or 2.3 a and 1.9 a confidence detections determined 
from Monte Carlo simulations). The combined significance of both lines is 3.0 tr, again determined 
from Monte Carlo simulations. The same observation yields a high statistics EPIC spectrum of 
the Coma cluster gas at the position of X Comae. We detect emission by Ne IX with a flux of 
2.5 ± 1.2 x 10~ 8 photons cm -2 s _1 arcmin~ 2 (90% confidence errors or 3.4 a confidence detection). 
These data permit a number of diagnostics to determine the properties of the material causing the 
absorption and producing the emission. Although a wide range of properties is permitted, values near 
the midpoint of the range are temperature ~4x 10 6 K, density ~ 6 x 10~ 6 cm~ 3 corresponding to an 
overdensity with respect to the mean of ~ 32, line of sight path length through it ~ 41 (Z/Zq)' 1 Mpc 
where (Z/Z & ) is the neon metallicity relative to solar. All of these properties are what has been 
predicted of the WHIM, so we conclude that we have detected the WHIM associated with the Coma 
cluster. 

Subject headings: galaxies: clusters: individual: Coma — intergalactic medium — quasars: absorption 
lines — large-scale structure of universe 



1. INTRODUCTION 

O ur current understanding of the state of cluster gas 
(e.g. iVoitlteOOl iBorgani et al.ll2005Tl requires some com- 
bination of preheating and radiative cooling of the am- 
bient gas, which is then further heated by the accre- 
tion shock as the gas falls into the cluster along with 
additional cooling and nongravitational heating after- 
wards. However our direct observational knowledge of 
the state of the gas before it is accreted onto the clus- 
ter is limited. Star formation, AGN activity and accre- 
tion shocks onto large-scale structures are all possible 
for the energy source of the preheating (see revi ews in 
e.g. IBorgani et al.ll2002l iDos Santos fc Dot& [2002). Fur- 
ther, the material that will later become the cluster gas is 
thought to be related to the low redshift missing baryons, 
most of whic h is suggested from recent numerical sim u- 
lations (e.g.. ICen fc Ostrikerl TT999t IDave et al.112001^ to 
reside in a warm-hot intergalactic medium (WHIM) with 
temperatures of 10 5 -10 7 K. Therefore, detecting warm- 
hot gas around clusters of galaxies is crucial to under- 
stand their formation as well as to settle the missing 
baryon problem. 

The soft X-ray excess above the harder intraclus- 
ter medium (ICM) emission reported for some clusters 
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(iLien et al.lll99ft iBowyer et al.lll999t iBonamente et alJ 
120021 iFinoguenov et alJ 120031: iKaastra et al.ll2003[) may 

be signaling a fortunate orientation of a filament 
containing WHIM at those clusters. In particular, 
IFinoguenov et ail l)2003f ) determined the WHIM density, 
temperature and abundance of heavy elements, assuming 
the soft excess in the EPIC spectra of the Coma cluster 
outskirts is due to a filament extending ~20 Mpc along 
the line of sight, corresponding to the excess galaxies in 
front of Coma. The conclusions are, however, strongly 
dependent on that assumption. 

An additional uncertainty common to all studies of low 
redshift cluster soft excess is the presence of emission 
lines in the Milky Way and even interplanetary fore- 
ground llMcCammon et alJ 120021 iWargelin et alJ 120041: 
iSnowden et alJl2004l) that are not easily separable from 
cluster emission given the spectral resolution of a CCD. 
So a confirmation of cosmological origin of the soft com- 
ponents in clusters is required. 

A direct way to confirm the existence of the warm-hot 
gas in cluster outskirts is to detect absorption lines in X- 
ray spectra of background quasars with a high resolution 
grating spectrometer. The spectral resolution of these 
instruments is sufficient to separate absorption by clus- 
ter WHIM from that of foreground contamination due 
to the interstellar medium in our Gala xy or the inter- 
planetary medium of our Solar System ijFutamoto et alJ 
120041 lYao & Wandl2005l) . Given the expected tempera- 
tures of this gas and the cosmic abundances of the ele- 
ments, the strongest lines should be the resonance lines 
of hydrogen-like and helium-like oxygen and neon. Addi- 
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TABLE 1 
Instrument mode 



Fig , 1. — An image of the C oma cluster in the 0.5-2.0 keV band 
from Finogucnov ct al. (2003). The coordinate grid is RA, Dec 
(J2000.0). 



tional information is available by combining absorption 
and emission measurements, particularly from the same 
ion. We can measure the density of the g well as 

constraining the geometry of the emitti ng zone directly 
ijKrolik & RavmondllT98alSarazinlll989fl . 

However desirable these measurement s are, they are 
at the limit of curren t instrumentation ijKravtsov et al.l 
l25o2lViel et all2003|) and detecting the WHIM will only 
be possible in special circumstances. The simulations 
mentioned previously indicate that the WHIM near clus- 
ters of galaxies should be denser than average. Therefore, 
we expect a higher absorption signal from it compared to 
that along random sight lines. One exam ple is a marginal 
detect ion of O VIII absorption line by iFuiimoto et al.l 
(2004) toward the Virgo cluster. It is possible to detect 
the average density WHIM in the X-ray spectra of ex- 
traordinarily bright background sources s uch as blazars 
parti c ularly when they are in outburs t. llNicastro et al. | 
2002t iFang et all l2002t iMathur et all 120031: iFang et al ' 
2003t IRasmussen et all I2003t IMcKernan et alJ I2004t 
NrcastToet^lTl2005f r However, even the most con 



vincing detection llNicastro et al .' 2005) is controversial 
l)Rasmussen et al.ll2006t iKaastra et alJ 20J36). The best 



case of a bright blazar behind a cluster is so far unknown. 
Consequently, there is still no firm X-ray detection of ab- 
sorption due to the WHIM. 

In this paper we present our XMM-Newton RGS obser- 
vations of the Seyfert 1 AGN X Comae, which is located 
behind the Coma cluster to the north but well inside the 
cluster virial radius. We also describe the EPIC spectra 
of the Coma cluster gas at the position of X Comae that 
were obtained simultaneously. We assume a Hubble con- 
stant of 70 km s _1 Mpc^ 1 or h^o = 1 and S7 m o = 0.3, 
^ao = 0.7. Unless otherwise stated, errors in the fig- 
ures are at the 68 % confidence level and at the 90 % 
confidence level in the text and tables. 

2. THE COMA CLUSTER, X COMAE AND OBSERVATION 
SUMMARY 



s ystem ic redshift is z coma = 
3dlll999(K which yields a scale of 
/lyp 1 Mpc per degree for our assumed cosmology. 



The Coma cluster 
0.0231 ijStruble fc Roo 
1.68 

The redshift dispersion of the cluster galaxies is a ga i 



Instrument 




Mode 


Filter 


RGS 




Spectro+Q 




EPIC pn 




Full frame 


Medium 


EPIC MOS 




Full frame 


Thin 






TABLE 2 




Exposure times and fluxes of X Comae 


Date 


Duration 


Net exposure 


Flux a (0.3-2.0 keV) 




ks 


ks a ks b 


ergs cm -2 s _1 


2004 Jun 6 


102.6 


71.5 60.1 


1.92 X 10~ 12 


2004 Jun 18 


108.3 


54.7 46.2 


1.48 X 10~ 12 


2004 Jul 12 


104.2 


45.6 39.5 


1.36 x 10~ 12 


2005 Jun 27 


55.9 


23.9 20.3 


1.02 x 10" 12 


2005 Jun 28 


80.8 


62.5 57.6 


1.46 x 10- 12 


Total 


451.8 


258.2 223.7 


1.66 x W)- rz 



a RGS 
b EPIC pn 



3.44 x 10~ 3 . T he integrated temp erature of the clus- 
ter is 8.21 keV ijHughes et all ^993). which implies the 
radius within which the cluster's average density is 200 
times the critical den sity, r?nn, is 2.4 hm Mpc. At a red- 
shift of 0.091±0.001 l|Branduardi-Ravmont et"all fl985^ . 
X Comae is the brightest X-ray object behind the Coma 
cluster. Its position is (13 h 00 m 22M7, +28°24'2".6) 
(J2000.0). As shown in Figure □ it is 28.4' or 0.79/1^ 
Mpc or 0.33r2oo north of the cluster center, defined to be 
at NGC 4874. 

We made five observations with XMM-Newton 
ijJansen et al.ll200l|) in 2004 and 2005 for a total of 451.8 
ks. The instrument mode and the filter used are shown 
in Tabled and the gross and net exposure times and flux 
of X Comae are summarized in Table El X Comae had 
a flux at or below its hitherto historical low during our 
observations. 

3. ABSORPTION LINES IN THE RGS SPECTRA 

We reduced the RGS <|den Herder et al.112001^ data us- 
ing the XMM-Newton Science Analysis System (SAS) 
version 6.1.0, with standard parameters. We checked the 
data from a source-free region of CCD 9 for background 
flares. The regions were CHIPX = 2 to 341, CHIPY = 
2 to 38 and 87 to 127 for RGS1 and CHIPX = 2 to 284, 
CHIPY = 2 to 37 and 87 to 127 for RGS2. We accu- 
mulated the signal photons only when the count rates 
in these regions were less than 0.1 s _1 (in PI range 80 - 
3000) for both RGS1 and 2. This rather severe threshold 
was determined empirically to yield the highest signal to 
noise background subtracted source spectrum. The net 
exposure time was about 60 % of the total exposure and 
is summarized in Table |21 

The spectra were extracted after merging the five data 
sets and then binned by a factor of four, resulting in a 
final bin width of 0.035 A at 11.5 A and 0.046 A at 23 A. 
These bin width are about half of the average FWHM 
wavelength resolution of the RGS (0.067 A). The red- 
shift widths corresponding to the resolution are 0.0058 
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Wavelength (Angstrom) Wavelength (Angstrom) 

Fig. 2. — (upper panel) RGS spectra around the Ne and O Ka lines. Black and red represent RGS1 and RGS2, while crosses (+) and 
triangles (A) show source plus background and background, respectively. The solid lines show the best fit model of source plus background 
and background. The wavelengths corresponding to redshifted Ne and O lines (z = 0.0231) are indicated with vertical dashed lines, (lower 
panel) Source plus background residuals from the best fit broken power law model. The model is continuum only; it has no absorption 
features. 



and 0.0029 at 11.5 A and 23 A respectively. The back- 
ground spectra were similarly produced from the same 
data sets using the pixels outside the region where source 
photons were dispersed. The background spectra were 
then scaled to account for the different areas used to ex- 
tract source and background photons. Since the number 
of photons is small, we used the C-sta tistic (maximum- 
likelih ood) method for model fitting. iRasmussen et all 
(2006) pointed out that merging different observations 
may introduce artificial absorption features in RGS spec- 
tra. However, their suggested features have about one 
order of magnitude smaller equivalent widths than those 
we discuss below. They are smaller than our statistical 
errors. We tried the procedure to search for bad columns 
that may cause artificial absorptions according to their 
Appendix C. It showed no apparent artificial features 
within our statistics. 



TABLE 3 

Results of fitting the RGS spectra of X comae with 

BROKEN POWER LAW MODELS 



Component 



Unit 



Value 



7V H 

Source T a (E < 0.75 keV) 
Source T a (E > 0.75 keV) 
Source Normalization 13 
Background T a (E < 0.75 keV) 
Background T a (E > 0.75 keV) 
Background Normalization 13 



9.3 X 10 19 (fixed) 
1 74+0.15 
A - ' -0.29 
2.41 ± 0.25 

(6.88±g;S) x 10- 4 



3.73Z 



U-statistic 
free parameters 
d.o.f. 



-0.66 
3+0.18 
-0.12 

1.96 ±0.23 

(0.79 ±0.08) x 10- 



708.63 
6 

628 



a photon index 

b in units of photons keV - 1 



cm 2 s 1 at 1 keV 



3.1. Detection of absorption features and their 
equivalent widths 

Figure shows the RGS spectra of X Comae in the 
Ne region (11.5-16.0 A) and the O region (18.0-22.7 A). 
We calculated the ratio of the data to a continuum only 
model in order to estimate the statistical significance of 
possible absorption features at the wavelengths of Ne IX, 
Ne X, O VII and O VIII Ka resonance lines of the Coma 
redshift (z = 0.0231) in a model independent way. That 
is, we calculated 

(source + background data) — (background data) 

Ratio = -. 

(source + background model) — (background model) 

(1) 

This procedure is called the ratio method in what fol- 
lows. 

The source model is a broken power law multiplied 
by the Galactic absorptio n of — 9.3 x 10 19 cm" 2 
l)Dickev fc Lockmanlll990|) . The background model is a 
different broken power law without Galactic absorption. 



In order to determine the continuum model, we fitted 
the source plus background and background spectra of 
RGS1 and 2 simultaneously using XSPEC version 11.3 
in the wavelength region 11.5 A-22.7 A. Regions around 
the Ne IX, Ne X, O VII and O VIII Ka lines, which are 
defined as z coma - 3 cr ga i < z < z coma + 3 <T gal , were 
excluded from the fit. 

The best fit model and parameters are shown in Fig- 
ure |2 and in Table respectively. The Ratios around 
Ne IX, Ne X, O VII and O VIII lines vs. z are shown in 
Figure □ The FWHM wavelength resolution of RGS is 
also shown as a horizontal line. We summed the RGS1 
and 2 data where both were available. Since the mapping 
from wavelength to redshift is a function of wavelength, 
we interpolated the counts to a common redshift bin size 
for all four lines. 

The three vertical dashed lines in Figure [3] indicate 
Zcoma and z coma ±2.5 x (jg a i. We defined redshifts within 
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TABLE 4 

Ratio for continuum and absorption spectral regions 





Continuum region a 


Absorption region a,b 


EW C 


Ne IX 


1.027 ±0.058 


0.782 ± 0.071 (98.0%) 


3.3 ± 1.8 eV 


NeX 


1.011 ±0.073 


0.950 ±0.092 (42.7%) 


0.8 (< 3.9) eV 


O VII 


0.908 ± 0.080 


0.927 ±0.103 (50.0%) 


0.7 (< 2.6) eV 


O VIII 


0.963 ± 0.054 


0.845 ± 0.071 (94.1%) 


1.7 ± 1.3 eV 


Avg. of Ne IX and O VIII 


0.993 ± 0.039 


0.813 ± 0.050 (99.7%) 





a Errors are quoted at 68% confidence level. 

b The probability that a simulated spectrum without absorption yields a smaller discrepancy from unity. 
c Errors are quoted at 90% confidence level, while upper limits are 2cr. 
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Fig. 3. — Ratios of the data to the continuum model (as defined 
in equation 1) versus z for Ne IX, Ne X, O VII and O VIII (top 
to bottom). Vertical dashed lines indicate z mma and ±2.5 X <r ga i. 
The average RGS1 and RGS2 instrumental resolution of 0.067 A 
(FWHM) corresponds to a redshift resolution of 0.0050, 0.0055, 
0.0031 and 0.0035 (top to bottom), which is indicated as a hori- 
zontal line at the lower right of each panel. 



z coma ± 2.5 x Cgai as absorption and the remainder in 
Figure as continuum, and then calculated the error- 
weighted average of the absorption and continuum ratios. 
The results are shown in Table 01 When the absorption 
redshift region was defined, we fixed its center to the 
apriori known z COIlm and chose its width to maximize the 
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Fig. 4. — Error- weighted average ratio of the data to the 
continuum model (as denned in equation 1) versus z for the 
Ne IX and O VIII lines. Vertical dashed lines indicate z CO ma 
and z C oma ± 2.5 X <T ga j. The two horizontal lines show the ±1 a 
error of the model continuum normalization with power law indices 
fixed at their best fit values. The average of the ratio in the six 
bins between the left and right vertical lines is 0.813 ± 0.050. The 
significance of this absorption is 99.7% according to Monte Carlo 
simulations. The average of the remaining ten continuum bins is 
0.993±0.039, which is consistent with unity. The FWHM averaged 
wavelength resolution of RGS is also shown as a horizontal line at 
lower right. 



Ne IX plus O VIII signal from among 2, 4, 6 or 8 binncd- 
by-four pixels; i.e., the region was determined after a 
four-trial optimization. 

Ne IX is the ion with the deepest absorption, and 
O VIII is the second deepest. The absorption ratio is 
below 1 for the Ne X and O VII lines as well, though 
they are not very significant. The continuum is always 
consistent with 1. This situation of strong Ne IX ab- 
sorption and weak absorption by the other three lines is 
often observed at much higher signal to noise in inter- 
stellar medi um features in the s pectra of galactic X-ray 
sources (e.g. lYao fc Wane| |2005). 

To improve the signal-to-noise, we made a grand error- 
weighted average of the ratios for Ne IX and O VIII lines 
and calculated the combined significance. The result is 
shown in Figure^lwhere the band around unity is +0.7%, 
— 1.0%, the la error of the model normalization with 
power law indices fixed at their best fit values. The aver- 
age of the grand-averaged ratio are also given in Table 0J 

Since the number of counts in each bin is not very 
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high (20-30 counts/bin), we investigated the significance 
of the absorption using Monte Carlo simulations. We 
made 1000 simulated spectra with no absorption which 
have the same statistics and the same response function 
as the actual data, and then calculated the ratios with 
the same procedure as above. That is, we calculted the 
ratio for Ne IX, Ne X, O VII and O VIII lines and grand 
error-weighted average of the ratios of Ne IX and O VIII. 
This calculation was done for 2, 4, 6 and 8 binned-by- 
four pixels around z coma , and then we chose the most 
significant grand average among the four trials. The sig- 
nificance of the absorption in our observed RGS data can 
be estimated as the probability that the simultated spec- 
tra without absorption yield a smaller discrepancy from 
unity. The significance for the grand error-weighted av- 
erage was 99.7%. We conclude that we have detected 
absorption by material associated with the Coma cluster 
with a significance of 99.7%. This significance is equiva- 
lent to 3.0cr of a Gaussian distribution. The significance 
of Ne IX and O VIII was 98.0% (2.3<r) and 94.1% (1.9a), 
respectively. 

The equivalent width EW of the absorption lines can 
be calculated from the ratios as (1 — Ratio) x AE, where 
AE is the energy corresponding to the width of 6 bins 
we used to calculate the ratios. The EWs for the four 
lines are also shown in Table Assuming that the 
absorption lines are not saturated, the column density 
A ion is calculated from EW as 



A io n = 



m e c(l 



EW 



Nu 



nhe 2 



9.11 x 10 15 cm 



fos 

,(l + z) EW 



(2) 



(3) 



fos leV 

l)Sarazinlll989t) where / os is the oscillator strength of the 
transition and the other symbols have their usual mean- 
ings. We used / os = 0.724 for the Ne IX, / os = 0.696 
for O VII, and f o s = 0.416 for Ne X and O VIII 
ijVerner et alJ 1*1996'). The column densities of the four 
lines are then iY Nc ix = 4.3 ± 2.3 x 10 16 cm" 2 , V NcX = 



16 



1.9 (< 8.8) x 10 

and iVovm = 3.7 ± 2.8 x 10 16 cm 



, AW 



0.9 (< 3.5) xlO 16 cm" 

,-2 



3.2. Fitting the absorption features and their equivalent 

widths 

Next we derived the equivalent widths of the four lines 
by another way: using model-fitting of the spectra. We 
adopted the same contimuum model as that in § 13.11 
and a boxcar profile to describe the absorption (NOTCH 
model in XSPEC). This absorption model multiplies the 
continuum by a factor of 



(1-F) 



for (A no tch — 

for all other, 



'•>$&■) < A < (A notch + ■ 



(4) 

where A no tch, W no tch an d F are the central wavelength, 
width and absorption factor, respectively. The equiva- 
lent width is given by EW — W notc h x F. We fixed A no tch 
to be the value corresponding the redshift of the Coma 
cluster. Since the significances of the absorption features 
were low, except for Ne IX, we fitted the Ne IX absorp- 
tion first and then fitted the other lines with Wnotch fixed 
to the best fit Ne IX value (in z, not wavelength). All to- 
gether the free parameters were: power law indices and 




0.005 0.01 0.015 0.02 0.025 
Width (converted to z) 



Fig. 5. — Contour plot of AC as a function of W no tch (con- 
verted to redshift) and F for the Ne IX line. Equivalent width is 
Wnotch x F - The f °ur contours are AC = 1.0, 2.71, 4.0 and 6.6 
corresponding to 68%, 90%, 95% and 99% confidence for one inter- 
esting parameter, respectively. These also correspond to la, 1.6a, 
2a and 2.6cr, respectively if the data were Gaussian distributed. 




0.01 



0.08 



0.03 



0.04 



Fig. 6. — Ratios of the data to the continuum model vsesus z 
for Ne IX (the same data as the top panel of Figure overlaid 
with boxcar absorption models. Red line shows the best-fit model, 
while green line displays (F, Wnotch) = (1-0, 3.1 X 10~ 3 ). Since the 
wavelength resolution of the RGS is comparable to the width of the 
absorptiion, the fit is not sensitive to F or Wnotch individually, but 
to the EW , which is proportional to F x Wnotch- Vertical dashed 
lines indicate 2 com and z„ a ± 2.5 x <r ga i. 



normalizations of the source and background continua, 
Wnotch and F for Ne IX absorption, and F for the other 
species. Data between 11.5-22.7 A were used in the fit. 

The best fit parameters are shown in Table U3 where 
Anotch and Wnotch are converted to redshift. The equiv- 
alent widths and their errors were determined taking 
into account the covariance between Wnotch and F, as 
is shown in Figure for Ne IX. Figure [5] indicates that 
this analysis is not sensitive to F or W no tch because the 
wavelength resolution of RGS is comparable to the width 
of the absorption. Figure HJ] compares the best-fit (red) 
boxcar absorption model and the model with F = 1.0 
and Wiotch = 3.1 x 10~ 3 (green). There is little differ- 
ence between them. Note that we can estimate the EW, 
the product of F and W no tch, more precisely than either 
F or Wnotch- From this figure, the equivalent width of 
Ne IX is estimated to be 3.7^2 ° e V (90% confidence er- 
rors). Best fit values and 2a upper limits are also given 
in Table for the Ne X, O VII and O VIII absorption 
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TABLE 5 

Results of fitting the RGS spectra of X Comae with broken power law plus boxcar absorption models 



Component 


Unit 


Value 


Continuum 


N H 

Source r (E < 0.75 keV) 
Source T (E > 0.75 keV) 


cm 1 


9.3 X 10 ly (fixed) 
-i 71+0.21 
1 -' J --0.26 
2.35 ±0.24 


Source Normalization 21 




6.92+°'^ x 10 -4 


Background r (E > 0.75 keV) 
Background T (E < 0.75 keV) 




3 64+ 11 
2.01±0.21 


Background Normalization 




(0.82 ±0.07) x 10~ 4 


NelX 13.447 A 


\ b 

''notch 

Wnotch 

F 




0.0231 (fixed) 
9.79 xlO -3 
0.41 


EW C 




3.7 ±5;2 eV 


iv ion 


cm z 


4.7 xlO 16 


NeX 12.134 A 


-^notch 

w notch h 

F 




0.0231 (fixed) 

9.79 XlO -3 (fixed to the NelX value) 
0.09 


EW C d 


eV 


0.8 (<4.7) 


N; c d 


_9 

cm 


1.9 (<10.5) xlO 16 


EW/EW Nclx d 

A r ion/A r NcIX d 




0.2 (<1.3) 
0.4 (<2.3) 


OVII 21.602 A 


-^notch 

w notch h 

F 




0.0231 (fixed) 

9.79 XlO -3 (fixed to the NelX value) 
0.04 


EW C d 


eV 


0.2 (<2.2) 


N; c d 

EW/EW Nolx d 

A r ion/A r NcIX d 


cm 


0.3 (<2.9) xlO 16 
0.06 (<0.59) 
0.06 (<0.63) 


OVIII 18.969 A 


\ b 

''notch 

Wnotch 
F 




0.0231 (fixed) 

9.79 XlO -3 (fixed to the NelX value) 
0.06 


EW C d 


eV 


0.4 (<2.2) 


N ion c d 

EW/EW KeIx d 

A r ion/A r NcIX d 


_9 

cm 


0.8 (<4.9) xlO 16 
0.10 (<0.59) 
0.2 (<1.0) 


C-statistic 




834.33 



Free parameters 11 
Degrees of freedom 745 

a In units of photons keV^ 1 cm~ 2 s _1 at 1 keV 
b Wavelength converted to redshift 
c Errors include covariance of W no t c h an d F 
d Upper limit is at 2 a confidence 



lines. The column density N- iaa of each ion, estimated 
from the EW are tabulated in Table as well. The de- 
rived equivalent widths are consistent within the errors 
with those obtained from the ratio of the data to the 
continuum model described in § 13.11 

The significance of the line was again calculated using 
Monte Carlo simulations. The C-statistic was improved 
by 7.57 when we added the Ne IX absorption line. The 
probability that the simulated data shows less improve- 
ment of the C-statistic was 99.2%, equivalent to 2.1a if 
the data were Gaussian distributed (compared to 98.0% 
from the ratio method). We took this value as the sig- 
nificance of Ne IX line with boxcar- fitting. In the case 
of boxcar-fitting, the O VIII line is less significant than 
that found with the ratio method (82.4% compared to 
94.1%) and hence combining Ne IX and O VIII absorp- 



tion did not improve the significance (98.5%) compared 
to that of Ne IX alone. The significances obtained with 
boxcar-fitting were slightly different from those by the 
ratio method, particularly for O VIII. This is probably 
because the two methods are not exactly identical: with 
boxcar-fitting the detector response is convolved with 
an assumed absorption shape (boxcar) that is same for 
the four lines, while no shape was assumed for the ratio 
method. 

4. EMISSION LINES IN THE EPIC SPECTRA 

We clean ed the EPIC X Comae d ata of flares as de- 
scribed by iFinoeuenov et al.1 120031). This proc edure 
yielded 223.7 ks of clean EPIC pn Intruder et alJl200lj) 
data. This exposure is among the longest ever made of 
a cluster with XMM-Newton. We considered the entire 
X Comae field, as well as dividing it into five concentric 
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Fig. 7. — ROSAT R4 band map (~0. 4-1.0 keV) near the Coma cluster jSnowden ct al. 1997) with analyzed regions of the EPIC overlaid. 
The R4 band map was smoothed with a two-dimensional Gaussian with a = 3 pixels. 



TABLE 6 

O AND NE LINE INTENSITY 



Field 



ObsID 



Exp. 
ks 



Distance from 
NGC 4874 
arcmin 



O VII 7 a 



O VIII / a 



NcIX I a 



X Comae 






26.5 


55.6±?;J 


15-4tM 


6.0 ±0.9 


X Comae 1 






17.6 




12.9^ 7 


11 5 +6 ' 4 


X Comae 2 






20.6 


59.7I&S 




7 7+2.8 
'•'-2.7 


X Comae 3 






25.9 


54.7-JI 




5.7 ± 1.7 


X Comae 4 






31.8 




16.0 ±2.1 




X Comae 5 






38.2 


53 


15.6^2.8 


5 2+ 1 ' 7 


Coma 


0124711501 


15.4 


46.2 


71.5 ±5.1 


28.1 ±3.4 


4.5 ± 1.8 


1253+275 


0058940701 


10.5 


71.6 


44.9jlg 7 


15.1 ±4.5 


1.1 (< 3.9) 


3C 284 


0021740201 


34.2 


156.4 


62.2 ±4.0 


26.4 ±2.5 


3.7 ± 1.2 


P Comae 


0148680101 


29.4 


162.7 


70.2 ±5.5 


17.5 ±3.2 


3.7 ± 1.6 


Average bkgd 






109.2 


62.2 ±9.8 


21.8 ±5.2 


3.3 ± 1.2 


Average bkgd b 






130.2 


59.1 ± 12.0 


19.7 ±5.5 


2.8 ± 1.4 


Wtd avg bkgd 






109.2 


63.6 ±2.5 


23.1 ± 1.6 


3.5 ±0.8 


Wtd avg bkgd b 






130.2 


61.1 ±2.9 


21.8 ±1.8 


3.3 ±0.9 



a Surface brightness in units of 10 8 photons cm 2 s 1 arcmin 2 

b Except Coma field, which may have a contribution from emission associated with the cluster. 



sectors of annuli, approximately centered on the center 
of the cluster. 

Determining a reliable background is crucial for this 
work because the temperatures of the WHIM and the 
Milky Way interstellar medium are similar and the soft 
emission in the X Comae field is not very bright. Many 
previous attempts to measure the soft em i ssion from the 
Coma cluster (e.g. lFinoguenov et alJl2003llKaastra et al.1 
l2003ll used ROSAT All-Sky Survey (RASS) data to ob- 
tain the values of the components of a galactic back- 
ground model. The resulting background subtraction 
is only as good as the model. Since then there have 
been several XMM-Newton observations serendipitously 
located around the Coma cluster, as we show in Figure[3 
Those observations provide a true background measured 
with the same instrument as our data and are therefore 
preferable for our analysis. We have examined the RASS 



data at the location of these fields to determine whether 
they are representative of the general background around 
the Coma cluster. One of them, HD111812 (ObsID 
0008220201), was located on an unusually bright spot of 
the RASS maps and has been omitted from the analysis. 
Some properties of the four remaining fields are shown 
in Table 

We performed the following analysis on these ten fields. 
We excised point sources including X Comae. We also 
excised diffuse sources using a wavelet-based detection 
algorithm, which were on a spatial scale of 4" to 4' and 
had a surface brightness comparable to that of the cluster 
emission. We extracted a spectrum from each of the ten 
fields and subtracted a filter wheel closed accumulation 
from each. We determined the vignetting correction as- 
suming a uniform source surface brightness in each field. 
We fitted the net vignetting corrected spectrum with a 
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Fig. 8. — EPIC pn spectrum from the entire X Comae field. 
Plotted is the ratio of the data to the smooth continuum with pa- 
rameters in Table |7| The grey line is a fit of three narrow width 
Gaussians to the residuals. The centers of the lower-energy Gaus- 
sians are fixed to O VII and O VIII at zero redshift and that of 
the higher-energy Gaussian is fixed to Ne IX at the Coma redshift. 
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Fig. 9. — Surface brightness of O VII and O VIII (crosses (+) 
and triangles (a), divided by ten) and Ne IX (circles (o)) versus 
distance from NGC 4874, which we take as the center of the Coma 
cluster. 



collisionally ionized thermal plasma model component 
(APEC in XSPEC) for the Coma hot gas, (except for 
the two background fields most distant from the Coma 
center), plus an APEC component for the Milky Way 
background, plus a power law component for the cos- 
mic X-ray background. The free and fixed parameters of 
the models are given in Table The energy range for 
the spectral fit was 0.4-7.0 keV, excluding the energies 
of strong detector lines (1.4-1.6 keV). We used the im- 
proved response matrix from the SAS v6.5 calibration re- 
lease. Our key finding is a clear detection of O VII, O VIII 
and Ne IX lines in excess of the continuum for the entire 
X Comae field, as shown in Figure |H1 Although there 
still remain small residuals at 0.97 keV, which might be 
identified as Fe XX L line complex, the statistical signif- 
icance of this feature is not high and hence we will not 
discuss it further. 
We then fitted narrow Gaussians at the energies of the 



O VII, O VIII and Ne IX lines. The results are in Figure 
121 an d Table [7| for the entire X Comae field and Fig- 
ure [5] and Table El for the other fields. The intensity of 
the Ne IX component increases towards the Coma clus- 
ter center (see Figure [5| and its intensity in the sector 
nearest the position of X Comae is above that of the 
background fields. The Coma intracluster medium is too 
hot to produce this line. We can also rule out a solar 
wind origin, as no variation with respect to the center of 
the Coma cluster in a single pointing is expected for that 
scenario. We thus conclude that we have detected Ne IX 
line emission from Coma cluster material that is cooler 
than most of the Coma intracluster medium. 

The intensity of the O VII and O VIII Gaussian com- 
ponents did not vary in an obvious way as a function of 
position (see Figure^ ■ The two oxygen lines show no en- 
hancement at the position of the cluster: their intensity 
in the five X Comae sectors is the same or even lower than 
that in the background fields. We thus conclude that all 
of the oxygen emission comes from the Milky Way soft 
background and not from material in the Coma cluster. 

The intensity of O VII and O VIII lines has a large 
scatter from one background field to another. The in- 
tensity in one of them (the 1253+275 field) is con- 
sistent with the soft X -ray background measured by 
IMcCammon et all (j2002D (they obtained 40.6 ± 10.9 and 
13.5 ± 5.4 xlO -8 ph cm -2 s -1 arcmin -2 , for O VII and 
O VIII, respectively), while the intensity in the other 
fields are larger than their values. On the other hand, 
the scatter of the Ne IX intensity among the back- 
ground fields is not as large and the average is consis- 
tent with the upper limit of IMcCammon et alJ l)2002f) . 
(< 5.4 x 10~ 8 ph cm~ 2 s _1 arcmin" 2 , from their Fig- 
ure 13). 

Do we expect not to be able to detect Coma cluster 
oxygen emission given the strength of the cluster neon 
line? The O VII surface brightness is 1.67 times that of 
Ne IX, for a temperature of 4 x 10 6 K and a Ne/O num- 
ber density ratio of 0.14 (see tj l5.1.1l for a justification of 
these values). This expected O VII line intensity in the 
entire X Comae field is approximately the dispersion of 
the other nine measurements (10.0 versus 8.0 xl0~ 8 ph 
cm~ 2 s arcmin -2 , respectively). Thus the different be- 
havior of the neon and oxygen emission lines is probably 
due to the much lower galactic neon background inten- 
sity that allows the Coma neon emission to be detected. 
The higher galactic oxygen background intensity masks 
the Coma oxygen emission in the X Comae fields. 

In the next section we will need the net Ne IX intensity 
at the position of X Comae. Of course this measurement 
requires an extrapolation to that position since the glare 
from the AGN prevents measurement of emission from 
the Coma cluster gas. We give in Tabled the numerical 
and error-weighted average of all four background fields 
and the three fields excluding Coma 0. The numerical 
average is more appropriate if there are real variations 
from field to field. Excluding Coma is appropriate if 
it has a higher intensity. Since none of the four averages 
are statistically distinguishable, we take for the Ne IX 
background the weighted average of all background fields 
since it has the lowest error. Similarly we take for the 
gross intensity the value from the entire X Comae field, 
since it is statistically indistinguishable from the sector 
closest to X Comae but has a lower error. The net inten- 
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TABLE 7 

Results of fitting the EPIC pn spectrum of diffuse gas at 

THE POSITION OF X COMAE a 



Component 


Unit 


Value 




Galactic absorption 


N H 


cm J 


8.0 X 10 iM (fixed) 


Coma hot 1CM 


kT 


keV 


, 7 r+U.3M 
• 3 - ,o -0.50 




Abundance 


solar 


0.47 ± 0.09 




z 




0.0231 (fixed) 




Normalization 13 




7.77±°:l° x 10 


-6 


Milky Way warm ISM 


kT 


keV 


u - 1 '^-0.001 




Abundance 


solar 


0, 1 (fixed) c 




z 




(fixed) 




Normalization 13 




5.5 ±0.2 x 10~ 


6 


Cosmic X-ray Background 


r 




1.40 (fixed) 




Normalization 11 




4.05±°; 2 g x 10 


-4 


O VII emission line 


E 


keV 


0.574 (fixed) 




r 




55.6^^ x 10" 


8 


O VIII emission line 


E 


keV 


0.654 (fixed) 




r 




15.4±^ x 10 - 


8 


Ne IX emission line 


E 




0.901 (fixed) 




I e 




6.0 ±0.9 x 10~ 


8 


Chi-squared 




1305.15 




Free parameters 




9 




Degrees of freedom 




1272 





a Entire X Comae field (456.9 arcmin 2 ) 

b J n c nndV/ATr(D\(l + z)) 2 per solid angle in units of 
10 14 cm -5 arcmin -2 , where D\ is the angular size distance to 
the source, n c is the electron density, njj is the hydrogen density 
and V is the volume. 

c Abundance of He, C, Fe, and Ni is fixed to 1.0, while that of all 
other elements, including O and Ne, is fixed to 0.0 

d In units of photons keV -1 cm -2 s -1 arcmin -2 at 1 keV 

e In units of photons cm -2 s — 1 arcmin -2 



sity is thus 2.5 ± 1.2 x 1CT 8 ph cm -2 s -1 arcmin -2 (90% 
confidence errors or 3.4cr detection). 

The intensity of Ne IX line is only 9% of the continuum 
level, as is seen in Figure |S1 Such a low intensity may 
mean that its measurement is subject to systematic er- 
rors in the continuum. We therefore repeated the entire 
preceding analysis with six additional continuum models 
that froze or thawed different components and/or added 
a soft proton component that was not folded through the 
mirror area. The dispersion in the seven measurement of 
the net Ne IX intensity was less than the above statistical 
error. We conclude that the statistical error accurately 
reflects the uncertainties of the measurement. 

5. DISCUSSION 

5.1. Properties of the warm-hot gas 

We have detected absorption features in the RGS spec- 
tra of X Comae at the redshift of the Coma cluster with 
a combined confidence of 99.7% (equivalent to 3. Oct of a 
Gaussian distribution) and line emission features in the 
EPIC pn spectra of diffuse gas at the position of X Comae 
at the 3.4ct confidence level. Although the significance of 
either one is not very high, the fact that we observed 
both absorption and emission from Ne IX at the Coma 
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Fig. 10. — The ratio of the theoretically expected EWs of Ne X, 
O VII or O VIII to Ne IX divided by the observed upper limits to 
that ratio. The dotted line shows NeX, while the dashed lines show 
the larger value for OVII or OVIII. Since the theoretical OVII or 
OVIII to NelX EW ratios depend on the number density ratio of 
Ne/O, we indicate three cases: 0.14 as the canonical value, 0.41 as 
the highest one found from literature, and 0.82 as the extremely 
high Ne/O case. The allowed temperature is the range in which 
both NelX and O curves are below 1. There is no allowed tempera- 
ture range for Ne/O = 0.14, while it is 2.0xl0 6 K < T < 5.8xl0 6 K 
for the extremely high Nc/O c ase. T he ionization fraction given in 
Table 2 and 3 of Chen ct al. (2003) for collisional ionization case 
was assumed. 



cluster redshift or position is additional support for the 
detection of this ion. In this section we give the proper- 
ties of the warm-hot gas that can be deduced from our 
observations under the assumption that the absorbing 
and emitting materials are the same and that material is 
uniformly distributed in a single phase in collisional ion- 
ization equilibrium. Although these assumptions have 
the virtue that they are simple and are consistent with 
the observations, it is entirely possible, even likely, that 
the actual situation is more complicated. In this case 
the properties we derive will be typical of the dominant 
phase of the material. 

5.1.1. Temperature and Ne/O number density ratio 

We can constrain the temperature of a plasma in ion- 
ization equilibrium from the ratios of absorption EWs 
between different ionization states of the same species. 
On the other hand, we can constrain the number density 
(or abundance) ratio of different elements from the ratio 
of absorption EWs between lines of different species. 

First, we investigate the allowed temperature and 
abundance for the boxcar-fitting case. Figure ^] plots 
the ratio of the theoretically expected EWs of Ne X, 
O VII or O VIII to Ne IX divided by the observed upper 
limits to that ratio. That is, it plots a ratio of ratios. The 
expected values were calculated using the ionization frac- 
tions of a plasma under collision al ion i zation equilibrium 
given in Figures 2 and 3 of lChen et alJ l)2003j) . In this cal- 
culation we assumed EW oc iVj on / os , i.e. no saturation. 
The dotted line shows Ne X, while the dashed lines show 
the larger value for O VII or O VIII. The allowed temper- 
ature is the range for which both the Ne IX and O curves 
are below 1. The upper limit of Ne X/Ne IX ratio gives a 
robust upper limit of the temperature, T < 5.8 x 10 6 K. 

The theoretical O VII or O VIII to Ne IX EW ratios 
depend on the number density ratio of Ne/O, while Ne/O 
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Fig. 11. — (left) Equivalent width ratio as a function of temperature T, based on Figure 2 and 3 of Chen ct al. (2003) in the collisional 
ionization equilibrium case. Solid and dashed curves correspond to OVII/OVIII and NeX/NelX, respectively. The horizontal lines represent 
the upper limits determined by this work. The allowed temperature is that for which the curved lines are below the horizontal upper limit 
lines of the same style, i.e., T > 2.2 X 10 6 K for OVII/OVIII and T < 5.7 X 10 6 K for NeX/NelX. (right) Ne/O number density ratio as 
a function of temperature T. The three curves correspond to the best-fit values and ±1<t errors, respectively, obtained from the observed 
NelX/O VIII ratio. 



ratio is not well known even for solar or interstellar val- 
ues. T he widely- used solar Ne/O number density ratio 
is 0.14 i|Anders k. Grevessdli989() . However it is difficult 
to measure the abundance of Ne due to its absence in 
the solar photospheric spectrum or in meteorites. Fur- 
ther, the O abundance also contains uncertainties; im- 
proved modeling of solar photosphere lines has yielded 
lowe r abundances for som e elements including O and 
Ne !|Asplund et all 12005^ . Recent measurements us- 
ing coronal X-ray or quiet Sun EUV lines find Ne/O 
< 0.30 (95% confidence) or 0.17±0.05 (estimated sys- 
tema tic errors) respectively l|Schmelz et al.ll2"005l lYound 
I2005D . A higher Ne/O ratio, 0.41, was suggested by 
X-ray sp ectroscopy of mostly giant stars and multiple 
systems (iDrake fc Testa 2005J). Recent measurements of 
sightlines to 4U1820-303 and LMC X-3 through the in- 
terstellar medium yield 0.20i ; 7 (9 0% confidence) and 
>0.14 (95% confid ence) respectively i|Yao fc Wanell2006t 
IWang et al.M2005|) . Thus the situation is not yet clear 
with reported values ranging from 0.14 to 0.41. 

Since our boxcar-fitting cannot determine the Ne/O 
ratio, we indicate three cases in Figure ITU1 0.14 as the 
canonical value, 0.41 as the highest one found in the 
literature, and 0.82 as an extremely high Ne/O case. 
We found that Ne/O > 0.25 (2a) is necessary in or- 
der to reproduce the boxcar-fit results with an assump- 
tion of single temperature plasma; for example, there 
is no intersection of O curve and 1 for Ne/O = 0.14. 
Thus, the canonical Ne/O ratio (Ne/O = 0.14) is re- 
jected in this analysis. On the other hand, if Ne/O was 
larger than 0.25, there was allowed temperature range. 
Even if adopting an extremely high Ne/O ratio, Ne/O 
= 0.82, the temperature range can be constrained to be 
T > 2.0 x 10 6 K. Combining this constraint with that 
from Ne X/Ne IX ratio, we obtained a conservative tem- 



perature range, 2.0 x 10 6 K < T < 5.8 x 10 6 K. 

We next used the EWs determined with the ratio 
method (S 13.11) . This method is independent of the in- 
trinsic shape of the absorption. The left panel of Fig- 
ure 1111 shows the theoretical equivalent width ratios of 
O VII/O VIII and Ne X/Ne IX as a functio n of temper- 
ature T again based on Figures 2 and 3 of iChen et alJ 
l|2003]l in the case of collisional ionization equilibrium. 
Note that here we use the O VII/O VIII ratio because 
O VIII was also detected at more than 90% confidence 
with the ratio method. We show the observed 2a upper 
limits obtained in § 13.11 with horizontal lines in the left 
panel of Figure 1111 The allowed temperature range is 
that for which the curved lines are below the horizon- 
tal upper limits of that ratio, i.e., T > 2.2 x 10 6 K for 
O VII/O VIII and T < 5.7 x 10 6 K for Ne IX/Ne X. Com- 
bining these two constraints, we restricted the tempera- 
ture of the plasma to be 2.2 x 10 6 K < T < 5.7 x 10 6 K 
at the 2a confidence level. 

From the Ne IX/O VIII ratio, we can investigate the 
allowed Ne/O number density ratio as a function of tem- 
perature T. This is shown in the right panel of Figure lTTl 
The ±1<7 error range is also shown. The allowed Ne/O 
ratio is not very sensitive to temperature; it is roughly 
0.3 ± 0.15 to 0.5 ± 0.3 within the allowed temperature 
range found above. Although our best-fit Ne/O ratio 
suggests a Ne overabundance, all of the above including 
the canonical value of 0.14 is within our ~ la confidence 
limit. Note that both temperature range and Ne/O ra- 
tio are consistent in the two estimates (boxcar fitting and 
ratio values). 

5.1.2. Density, line of sight length and metallicity 

We can constrain the average hydrogen density nn of 
the warm-hot gas and its path length L along the line of 
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Fig. 12. — Derived ripj (upper) and ZL (lower) as a function of 
T. Solid and dashed curves represents best-fit values and ler con- 
fidence regions, respectively. Vertical dotted lines indicate all owed 
temperature range, 2.0 X 10 6 K < T < 5.8 X 10 6 K (see FigurclTTl. 



sight by combining absorption and emission observations 
of the same ion, Ne IX in our case. The column density 
of an ion is 

N ion = fionZnnL (5) 

where Z is the abundance of the element relative to hy- 
drogen. Equation (2) gives N ion as a function of the ob- 
served EW . The surface brightness of an emission line / 
is 



C 



(1 



\3 H 



(6) 



where the exponent of (1 + z) is three instead of four 
because we measure surface brightness in photons, not 
ergs, and C is a coefficient depending on temperature 
but not on the elemental abundance. Solving these two 
equations simultaneously gives 



n H 



/ion(l + z) 3 I 



and 



ZL = 



C 



c 



N 2 



/L(i + *) 3 i 



The coefficient C is 

1 / n c 
4ir \riji 



C 



1 V j? 
Z{M) V E i ' 



(7) 
(8) 

(9) 



Here n e /nn = 1.1 7, Z(M) is the abu ndance of the el- 
ement assumed by iMewe et all (119851) (8 32 x IP" 5 fo r 
neon), P' is tabulated in Table IV of lMewe et alJ (|1985[) . 
E is the energy of the line and the sum is over all 
lines not resolved with CCD energy resolution. For the 
Ne IX emission C comes from six lines: Ne IX res- 
onance line (13.44 A), three satellite lines of Ne VIII 



(13.44 A, 13.46 A and 13.55 A), Ne IX intercombination 
line (13.55 A) and Ne IX forbidden line (13.70 A). 

We derive the parameters of the material producing 
the Ne IX absorption and emission as a function of tem- 
perature from its column density and intensity at the 
position of X Comae found at the ends of Sections 13.21 
and 01 respectively. For example, a temperature of 



4.0 x 10 6 K gives C = 5.37 x 10" 
and / ion = 0.50. Then 
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photons cm s sr 



(Z/Z e )L = 41 Mpc 



2.5 x 10- 8 J V 4.4 x 10 16 cm- 2 
/ \ " 1 ( JVNcIX 



2.5 X 10-8 



4.4 X 10 16 cm- 2 
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where Zq = 1.23 x 10~ 4 is the solar abundance of Ne 
relative to H l)Anders fc Grevessell989D . The derived val- 
ues of rifj and ZL depend strongly on the temperature. 
We show in Figure El the values as a function of tem- 
perature, where solid and dashed curves represent best- 
fit values and la confidence regions, respectively. Verti- 
cal dotted lines indicate the allowed temperature range, 
2.0 x 10 6 K < T < 5.8 x 10 6 K (see Figure E3J. We con- 
strained riH to be 2 x 10~ 6 cm~ 3 < tih < 8 x 10~ 5 cm~ 3 
and ZL to be 1 Z Mpc < ZL < 300 Z Q Mpc. These 
limits are at approximately the 3cr confidence level since 
they sum the 2a limit on temperature and the la limit 
on riH and ZL (i.e., do not sum the errors in quadra- 
ture). The derived hydrogen density corresponds to 
10 < 5 < 400. Here S is the overdensity with respect 
to the mean density of the universe: S = nn/n~}i with 
Uh = XQ,bp c (l + z) 3 /mn- X = 0.71 is the hydrogen- 
to-baryon mass ratio and fib = 0.0224ft, -2 is the baryon 
density relative to the critical density p c . 

Although the constraint we obtained on ZL is 
not tight, we note that the lower limit on L is > 
1 Mpc {Z/Zq)- 1 . The measured value of Z/Z Q at the 
position of X Comae is 0.16±0.05 (6 8% confidence error) 
for iron Ipe Grandi fc Molendil200H) where Z fl) = 4.68 x 
10~ 5 is the solar abundance of iron l)Anders fc Crevessel 
1989). Thus, the scale of the plasma that contains 
Ne IX is at least 2.6 times the size of the Coma cluster 
(?*200 = 2.4/itq 1 Mpc) if the Ne and Fe abundances are 
approximately the same. For at least two clusters (2A 
0335+096 and Sersic 159-03) RGS observations show 
that the Ne/F e is a pproximately unity or slightly less 
l|Werner et alJl200fit IdePlaa et al.l l2006). We conclude 
that the material containing Ne IX is likely not in the 
Coma cluster. 

5.2. Can the Ne IX emission come from material within 
the cluster virial radius? 



We investigate the suggestion bv I Cheng et alJ lj2005T) 
that gas associated with merging sub-groups inside the 
cluster virial region, which preserves its identity for a 
while before being destroyed by the hot intracluster 
medium, is responsible for the cluster soft excess. In par- 
ticular, we determine whether this material can produce 
the Ne IX emission we observe. We disregard the con- 
straints from the absorption measurements here. Rather 
we assume that the temperature of this material is 2 x 10 6 
K, the peak of the Ne IX ion fraction, and that the ma- 
terial is in pressure equilibrium with the ICM. Neither 
of these assumptions is very constraining. Changing the 
temperature to 4 x 10 6 K, the midpoint of the allowed 
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range found previously, changes the results calculated be- 
low by less than a factor of two. The second assumption 
yields a density of warm-hot material similar to the dens- 
est regions of groups, which we could have assumed at the 
outset since that is the suggestion we are investigating. 

First we need the properties of the hot ICM. The 
emission- weighted temperature at the posi tion of X Co- 
mae is 7.4 keV from the temperature map of lHonda et al.l 
(1996). The emission measure- weighted hydrogen den- 

sity at the positi on of X Comae is 2 x io- 4 /i^ 2 cm 
ijBriel et al.l 11992^) . Next we need the properties of the 
warm- hot material. A temperature of 0.172 keV and 
pressure equilibrium with the ICM yields a hydrogen 

density of 8.6 x lO^h 1 ^ 2 cm" 3 , or 5 ~ 3.4 x 10 4 , and 

entropy ~ 4 keV h 7 ^ 3 cm 2 . The latter two values 
place the warm-hot material on the extre me high den- 
sity ta il of the phase plot in Figure 4 of iCheng et al.l 
(2005). Equation (8) with C(2 x 10 6 K) = 1.19 x 
10 -13 photons cm 3 s _1 sr _1 , the above hydrogen density 
and the observed Ne IX surface brightness imply the path 
length through the material L is ~ 100 pc. The sound 
crossing time across L is ~ 5 x 10 5 years, during which 
the material travels ~ 600 pc moving at a ga i c. Thus 
pressure equilibrium is a good assumption because the 
ICM pressure hardly changes over such a small distance. 

How long can this warm-hot material survive? We as- 
sume its size across the line of sight is also L, that is 
the warm- hot material is tiny blobs of diameter ~ 100 
pc. Following iCowie & McKeei l)1977[) . the classical heat 
conduction across the blob-cluster interface is saturated. 
The saturated evaporation time is ~ 1 x 10 4 years, from 
their equation (64). Any blobs that might exist are very 
quickly destroyed. Further they are not replenished since 
the time for group gas at a temperature of 2 keV and the 
above density to cool to a temperature at which there is 
a significant population of Ne IX is ~ 5 x 10 10 years. 

Even after the blobs evaporate it takes some additional 
time for the Ne ion distribution to equilibrate to that ap- 
propriate to the new temperature in which the Ne finds 
itself. The longest lived ion capable of emitting the Ne IX 
resonace line we detect is Ne X, which it does by electron 
capture to an excited level followed by radiative deexcita- 
tion. The ionization time Ti on to convert Ne X to Ne XI is 
mostly determined by the collisional ionization efficiency 
<SNex; i.e., Ti on ~ l/^e^NeX- The recombination rate is 
negligibly small and the ionization time to convert Ne IX 
to Ne X is smaller than that for Ne X to Ne XI in our 
case. The empirical formula f or SueX f° r coronal plasma 
is given m Mc Whirteil l)1965l equation 40 of Chapter 5) 
as 

S NoX = 1.10 X lCr 5 J^ Tc/x)1/2 exp (--*-) cm' 3 S ~\ 
X 3/2 (6 + kT c / x ) V kTj 

(12) 

where T and \ are the electron temperature and the ion- 



ization energy in eV, respectively. Substituting the prop- 
erties of the hot gas, T c = 7400 eV, \ = 1360 eV, and 
n c = 2.3 x 10~ 4 h^cm -3 , we find Ti on = 3.7 x 10 6 years. 
The evaporation and the equilibration times are much 
smaller than the cluster crossing time, the characteristic 
scale of the situation. Thus we conclude that, under rea- 
sonable assumptions, material within the cluster virial 
radius is not capable of producing the observed Ne IX 
emission unless we are viewing the Coma cluster at a 
very special epoch. 

6. SUMMARY 

Our main result is the detection of absorption and 
emission lines from Ne IX associated with the Coma clus- 
ter of galaxies. The absorption is centered on the previ- 
ously known Coma redshift, although it is of marginal 
statistical significance (98.0% to 99.2% confidence de- 
pending on the analysis method). The emission is sta- 
tistically significant (3.4cr) and is positionally coincident 
with Coma. While the absorption line is resolved, its 
width is much smaller than the spectral resolution of 
the emission line data. We do not know if the emission 
and absorption lines have the same profile. The prop- 
erties of material causing these lines, assuming it is in 
a single phase in collisioinal ionization equilibrium, are 
constrained as follows: temperature T is 2.0 x 10 6 K < 
T < 5.8 x 10 6 K, density n H is 2 x 10" 6 cm~ 3 < n H < 
8 x 10 -5 cm -3 , overdensity 8 is 10 < S < 400, line of sight 
path length through it is L > 6 Mpc (Z/0.16 Z Q )- 1 . 

These properties are similar to those expected of the 
warm-hot intergalactic medium and we conclude that we 
have detected it. The warm- hot intergalactic medium 
is also expected to be distributed in filaments connect- 
ing clusters of galaxies. Since X Comae lies behind the 
apparent continuation of the Coma-A1367 chain of galax- 
ies to the northeast of the Coma cluster l|Gavazzi et alJ 
1999), the material we detect may reside in this previ- 
ously identified filament. 
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